K-essential Phantom Energy: Doomsday around the Corner? 
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Abstract 

We generalize some of those results reported by Gonzalez-Diaz by further tuning the 
parameter (/3) which is closely related to the canonical kinetic term in fc-essence formalism. 
The scale factor a(t) could be negative and decreasing within a specific range of j3 (= —1/u), 
lo : the equation-of-state parameter) during the initial evolutional period. 
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1 Introduction 

The standard (hot) big bang (SBB) theory is an extremely successful one, and has been around 
for over 60 years, since Gamow originally proposed it [1-2]. Remarkably, for such a simple idea, 
it provides us with an understanding of many of the basic features of our Universe. All that 
you require in the cooking pot, are initial conditions of an expanding scale factor, gravity, plus 
the standard particle physics we are used to, to provide the matter in the Universe. However, 
as mentioned the hot big bang theory can successfully proceed only if the initial conditions are 
very carefully chosen, and even then it only really works at temperatures low enough, so that 
the underlying physics can be well understood. The very early Universe is out of bounds, yet 
there is a hope that accurate observations of the present state of the Universe may highlight 
the types of process occurring during these early stages, providing an insight on the nature of 
physical laws at energies which it would be inconceivable to explore by other means. Another 
unresolved issue is the cause of the apparent acceleration of the Universe today, as seen through 
the distribution of distant Type la Supernovae. 

The physical properties of vacuum phantom energy are rather weird, as they include violation 
of the dominant-energy condition, P + p < (with the equation of state : P = ujp [3-4], uj < — 1; 
P is the homogeneous dark energy pressure and p is the energy density), naive superluminal 
sound speed and increasing vacuum-energy density. The latter property ultimately leads to 
the emergence of a singularity usually referred to as big rip in a finite time in future where 
both the scale factor and the vacuum-energy density blow up [5-7]. For instance, Gonzalez-Diaz 
showed that a cosmological model with a singular big rip at an arbitrary finite time in the 
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future can be also obtained when the scalar field satisfies equivalent phantom-energy conditions 
in the case that it is equipped with non-canonical kinetic energy [4] for models restricted by a 
Lagrangian of the form C = K(cf))q(y) [5]. He claimed that, one can play with the arbitrary 
values of the prefactor ao for the scale factor expression in Eq. (10) and those unboundedly 
small positive values of t which satisfy the observational constraint [5] (note that in this model 
(3 = —l/uj^ = — 1/cj) and the currently observed cosmic acceleration rate [8] to check that such 
a set of present observations [8] is compatible with unboundedly small positive values of t (the 
assumption is, if a quintessential scalar field 4> with constant equation of state P§ = to^p^ = up^ 
is considered, then phantom energy can be introduced by allowing violation of dominant energy 
condition, P^ + p^ < 0). Note that, in general k-essence is defined as a scalar field with non- 
canonical kinetic energy, but usually the models are restricted to the above Lagrangian form 
(C) [9]. Normally we use units where 8nG/3 = 1 [4]. 

Gonzalez-Diaz proposed that, from the perfect-fluid analogy, we have for the pressure and energy 
density of a generic /c-essence scalar field </> [9] 

P^y) = K(cf>)g(y)/y, p^y) = -K (</>)</ (y), 

where 

g(y) = By?, 

with B and (3 being given constants such that B > and < (3 < 1 and the prime means 
derivative with respect to y. Thus, the equation-of-state parameter reads 

yg'(y) ' 

We set next the general form of the function g(y) when we consider a phantom-energy fc-essence 
field; i.e., when we introduce the following two phantom energy conditions: K((f)) < and 
P<t>(y) + P<t>(y) < 0- The weak energy condition p${y) > is presumed to be valid here so that 
g'(y) > 0. From the second phantom-energy condition, we then deduce that g(y) > yg'(y) (and 
g(y) > 0). Therefore the function g(y) should be an increasing concave function, that is we must 
also set d 2 g(y)/dy 2 < 0. 

Following [5], we then specialize in the case of a spatially flat Friedmann- Robertson- Walker 
spacetime with line element ds 2 = —dt 2 + a(t) 2 dr 2 , in which a(t) is the scale factor. In the case 
of a universe dominated by a fc-essence phantom vacuum energy, the Einstein field equations are 
then [4-5] 

3H 2 =p4y), 2H + P4> (y) + P^y) = 0, 

with H = a/a, the overhead dot meaning time derivative. 

By combining the two expressions above and using the equation of state he can obtain for the 
function g{y) as given above 
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and then the solutions (for his spatially flat case) 

a ^ = { t- U )%nw-w 0< ^^ (2) 

where i* is an arbitrary and positive constant (this positive i* solution family does represent an 
accelerating universe if a(t) increases once t increases [5]). Gonzalez-Diaz presented the behavior 
of a(t) in Fig. 1 (II) therein [5] and made remarks : of quite greater interest is the choice i* > 
for which the universe (Fig. 1(H), solid line [5]) will first expand to reach a big rip singularity 
at the arbitrary time t = t* in the future, to thereafter steadily collapse to zero at infinity; that is 
it matches the behaviour expected for current quintessence models with to < — 1. The potentially 
dramatic difference is that whereas in quintessence models the time at which the big rip will occur 
depends nearly inversely on the absolute value of the state equation parameter, in the present 
k-essence model the time t* is a rather arbitrary parameter [5]. 

The present author, however, based on similar derivations (1/y 2 = <ft 2 /2 [4-5]) could obtain Eq. 
(10) the same as that in [5] but different behavior of a(t) for t < region with (3 = 1/2223, 
0.60 and 297/299 (compared to the (II) of Fig.l in [5] which is for < (3 < 1). Our results are 
illustrated in Fig. 1 with (3 = 1/2223, 0.60, 0.85, 297/299, respectively (t* = 2.0). We can clearly 
observe that the curves of (3 = 1/2223,0.60,297/299 (in t < t* region) are different from that 
claimed in Fig. 1 (II) in [5]. Under these situations, values of a(t) (which are negative, if we 
define 

Bp 3(1- (3) 3(1 + 1) 3( W + 1)' {6) 



we have B p = 1/3333,1,99 for (3 = 1/2223,0.60,297/299, respectively; then (the scale factor) 



should be negative for t < tj) firstly decrease as (time) t increases until t — > i*. Here, ao is an 
arbitrary integration constant. 

To be specific, considering (3 = 0.60 and 297/299 (or co = -5/3 and -299/297), of which the 
role of dark energy can be played by physical fields with positive energy and negative pressure 
which violates the strong energy condition p + 3P > 0(lo > —1/3), we thus obtain 

a(t) = a (t-t*)-\ a(t) = _ ° (5) 

and there is no doubt that a(t) should be negative once t < i*. These cases correspond to the 
earlier collapsing ones (but not collapse to zero) instead of expansion. 

The (3 = 0.85 curve (in t < t* region), however, resembles that presented by Gonzalez-Diaz in 
Fig. 1 (II) of [5]. It means, for some cases of (3, the universe will not first expand to reach 
a big rip singularity at the arbitrary time t = t* in the future as Gonzalez-Diaz claimed in [5]. 
Cases of a(t) being negative (in t < t* region) should be interpreted in another way (at lease 
for cases of (3 = 1/2223,0.60,297/299 where the power : B p is, al least, an odd integer or the 
denominator of B p is an odd integer while the numerator of B p is normalized to be 1.)! 
In fact, we can set 
n 

B p = —, n,m^0, (6) 

m 
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where n,m are positive integers. We already demonstrated cases of B p = 1/3333,1, and 99 
which correspond to (n,m) = (1, 3333), (1, 1), (99, 1) or (3 = 1/2223,3/5,297/299 (or u = 
—2223, —5/3, —299/297). To let B p be a positive odd integer, we can select m = 1 and B, p = n. 
It leads to 

= 3n = Wp = — (7) 
3n + 2 2 + 3B p co ' { ' 

Under this choice, a(t) will decrease (for t < t*) following 

a(t) = (t=t^ = (^F (8) 

since n is an odd integer (t — i* < 0)! For example, with m = 1, n = 101, we have decreasing 
o(t) for /3 = 303/305! 

On the other hand, once we choose n = 1, together with m is a positive integer, it then sets 
Bp = 1/m and 

With this, as we have shown above (n = l,m = 3), a(i) will decrease for t < t* following 
/ 1\ _ Q o a _ a 

One another example, for (3 ~ 0, is m = 1001 (i? p = 1/1001) or /? = 3/2005, which leads to 

a 



a(t) 



To make a brief summary, for those cases considered in Eqs. (7) and (9), we have earlier 
collapsing cases (a(t) decreasing once t < t*; a negative divergence!). However, as we have 
illustrated above, cases which don't belong to the restriction of Eqs. (7) or (9), say, n is an 
even integer (m is still an odd integer), e.g., = 0.85 = 17/20 gives n = 34 and m = 9 or 
B p = 34/9, will match the results claimed in [5] : the universe will first expand to reach a big 
rip singularity at the arbitrary time t = t*\ At least, our results can also be extended to the 
cases considered in [3]. Acknowledgements. The author is partially supported by the Starting 
Funds of XJU-2005-Scholars. 
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P= 1/2223, co= -2223 

-€)■■ P=0.60, oo= -5/3 
-x- P=0.85, co= -20/17 
-V- p=297/299, co= -299/297 




Fig. 1 Evolution of the scale factor a(t) with cosmological time t for a function g(y) 
with the form given as By@ . If we choose that t* to be positive, then the constant t* (=2.0) 
becomes an arbitrary time in the future at which the big rip takes place. All units in 
the plot are also arbitrary. Cases of (3 = 1/2223,0.60,297/299 are different from that in 
Fig. 1 (II) in [5] therein. a(t) < in t < t* region and a(t) decreases as t increases 
until t — ► £*. The equation-of-state parameter to = —1/(3. 



